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This article is concerned with statistics of addition spectra of certain 

many body systems of identical particles. In the first part, the pertinent sys- 
tem consists of N identical particles distributed among K < N independent 
sub-systems, such that the energy of each sub-system is a quadratic function 
of the number of particles residing on it with random coefficients. On a large 
scale, the ground state energy E{N) of the whole system grows quadratically 
with N, but in general there is no simple relation such as = aN + bN'^. 
The deviation of E{N) from exact quadratic behavior implies that its second 
difference (the inverse compressibility) xn = E{N -|- 1) — 2E{N) + E{N — 1) 
is a fluctuating quantity. Regarding the numbers XN as values assumed by a 
certain random variable Xi obtain a closed- form expression for its distri- 
bution F[x)- Rs main feature is that the corresponding density P{x) = 
has a maximum at the point x = 0. As — > oo the density is Poissonian, 
namely, P{x) e~^. 

This result serves as a starting point for the second part in which cou- 
pling between subsystems is included. More generally, a classical model is 
suggested in order to study fluctuations of Coulomb blockade peak spacings 
in large two-dimensional semiconductor quantum dots. R is based on the 
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electrostatics of several electron islands among which there are random in- 
ductive and capacitive couplings. Each island can accommodate electrons on 
quantum orbitals whose energy depend also on an external magnetic field. 
In contrast with a single island quantum dot where the spacing distribution 
between conductance peaks is close to Gaussian, here the distribution has a 
peak at small spacing value. The fluctuations are mainly due to charging 
effects. The model can explain the occasional occurrence of couples or even 
triples of closely spaced Coulomb blockade peaks, as well as the qualitative 
behavior of peak positions with the applied magnetic field. 
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I. STATISTICS OF ADDITION SPECTRA OF INDEPENDENT QUANTUM 

SYSTEMS 



A. Motivation 

Statistics of spectra is an efficient tool for elucidating properties of various physical 
systems. So far, most of the effort is focused on the study of energy levels of a system with 
a fixed number of particles. In this context, one of the central earlier results is that the 
spectral statistics of many-body systems such as complex nuclei agree with the predictions 
of random matrix theory On the other extreme, it was found that level statistics of a 
single particle in a chaotic or disordered system also obeys a Wigner-Dyson statistics @,|]. 

Recently, experiments are designed to get information on the statistics of the addition 
spectra of electrons in quantum dots |^. The pertinent energy levels E{N) are the ground 
state energies of a system consisting of N electrons residing on a quantum dot, which is 
coupled capacitively to its environment. 

Let us single out two properties of the addition spectra of quantum dots. The ffist one 
is that, on a large scale, the energy E{N) grows quadratically with A^, while the second one 
is a consequence of charge quantization, namely, there is, in general, no simple relation such 
as E{N) = aN + bN'^. In this context, an appropriate quantity whose statistics is of interest 
is then the inverse compressibility, 

XN = E{N + 1)-2E{N) + E{N-1). (1) 

It is the deviation of E{N) from exact quadratic behavior which makes its second difference 
Xn non-constant. Indeed, in a recent experiment on large quantum dot ^ it was found that 
the inverse compressibility vanishes for numerous values of electron number N. 

In the present section we study the statistics of the addition spectrum of a simple physical 
system with the two basic properties mentioned above. One example of such a system 
is motivated by considering the electrostatic energy of large quantum dots (although it 
should be mentioned that the model is too simple to describe the actual physics, see next 
section). To be specific, we have in mind a system of K metallic grains such that the 
number of electrons on the i*'' grain is rij [i = 0,1,2, K ~ 1) and their sum equals N. 
The electrostatic energy of the pertinent system is a bilinear form in the numbers with 
a K X K matrix w = \C~^. Here C is a positive-definite symmetric matrix of capacitance 
and inductance coefficients. If the metallic grains are very far apart, the matrix C is nearly 
diagonal. Thus, we concentrate on the special case C = diag[Cj], for which the energy of 
the system is given by 
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K-l ^ 

E{N) = mill ^ —nl (subject to Y.f=t ni = N). (2) 

The minimum in (|]) is taken over all possible partitions of A^. 

Another example is the energy of a system composed of K different harmonic oscillators, 
among which one distributes spinless fermions. If there are fermions on oscillator 
i (whose frequency is Ui), then the energy of this oscillator (up to a constant) is Ei = 
huiTiiirii + 1), and hence the ground state energy of the system is 

K-l 

E{N) = min ^ E,, (subject to E^o' ni = N). (3) 

i=0 

Discussion is concentrated on the first example, which is borrowed from the electrostatics of 
quantum dots (^, and refer to the constants Ci as capacitors. A remark pertaining to the 
second example (the system of oscillators (^) are also presented. 

Regarding the numbers xn of ([1|) as values assumed by a certain random variable, the 
distribution of this random variable is the main focus of the present work, which culminates 
in Theorem 1, where we find a closed-form expression for the distribution. 

The problem of elucidating the (addition) spectral statistics of a a many-body system, 
consisting of several independent sub-systems (whose dependence of E on is known), 
looks deceptively simple. As will be evident shortly, this is not the case, and finding the 
distribution in question is quite a non-trivial task. Note that, even for a single particle 
system composed of several independent sub-systems {e.g., a system of a particle in several 



boxes), the derivation of level statistics requires a large degree of mathematical effort pi 



The rest of this section is therefore devoted to a rigorous derivation of our main results. 
We offer our apology to the causal reader who might find this part rather mathematically 
oriented. He or she is adviced to start looking at problem 1 below and leave this section 
after reading the statement of Theorem 1. 



B. Formalism 

Definition 1. Let {6n)'^=i be a sequence of real numbers and F a distribution function. 
The sequence {6n) is asymptotically F- distributed if 

|{l<r2<M:^^„<a;}| ^, , 

T7 ^ -T X 

for every continuity point x of F (where IS"! denotes the cardinality of a finite set S). 

An equivalent condition is the following. Denote by 6t the point mass at t, and let fi be 
the probability measure corresponding to the distribution F (namely, fi{A) = J lA_dF{x) for 
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any Borel set A). Then is asymptotically F-distributed if 

1^ (So, + 6e, + . . . + 6eJ ^j^^fi 

(the convergence being in the weak*-topology). 

The notion of asymptotic distribution has a stronger version whereby, instead of requiring 
only that initial pieces of the sequence behave in a certain way, we require this to happen 
for any large finite portion of the sequence. This leads to 

Definition 2. In the setup of Definition 1, is asymptotically well F- distributed if 

\{L<n<M ■.en<x}\ ^ ^ 

for every continuity point x of F. 

Recall that the density of a set A C N is given by 

D(A) = lim \^nhm 

^ ' M^oo M 

if the limits exists. If, moreover, the limit 

BDM = lim J ^— — ^ 

^ ^ M-L^oo M - L 

exists (in which case it is certainly the same as D{A)), then it is called the Banach density 
of A(cf. ip.72]). 

The following lemma is routine. 

Lemma 1. Let {9n)'^=i be a sequence of real numbers. Suppose N = Uj=i^i! where the 
union is disjoint. Let (6*^^^)^^^ be the subsequence of [On], consisting of those elements 9^ 
with n & Aj, 1 < j < r. 

1. If each (9^^) is asymptotically Fj -distributed for some distribution functions Fj, 1 < 
j < r, and T){Aj) = dj, 1 < j < r, then {9n) is asymptotically F -distributed, where 
F = Y.Ud,F,. 

2. If each (9!^^) is asymptotically well Fj -distributed and BD{Aj) = dj, then {9n) is 
asymptotically well F -distributed. 

Obviously, a general sequence on the line does not have to be asymptotically distributed 
according to some distribution function, but one would expect it of sufficiently "regular" 
bounded sequences. In our case, one might expect xn to be distributed according to some 
distribution function corresponding to a measure centered at about 1/C. However, this is 
not the case. In fact, the measure in question is supported on a finite interval, and is a 
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convex combination of an absolutely continuous measure with decreasing density function 
on some interval [0, a] and the point mass 5a at the right end a of that interval. 
We have defined E{N) indirectly by means of the following 

Problem 1. For each non-negative integer N, find non-negative integers 
no, rii, . . . , riK-i, satisfying no + rii + . . . + nx-i — N, for which Z^^q^ ' minimal. 

It turns out that this problem is intimately related to a second optimization problem. 
Put Wi — -^,, < i < K — 1, and let A denote the set of all positive odd multiples of the 
numbers 

A = [wq, 3-^0, Swo, . . . , wi, 3wi, Swi, . . . , wk-i, "iwK-i, 5wk-i, • • •} • 

Here we treat A as a multi-set, or a sequence, in the sense that if some elements appear 
in this representation of A more than once (which occurs iff some ratio Wi/wj is a rational 
number with odd numerator and denominator) , then we consider A as having several copies 
of these numbers. 

Problem 2. For each non-negative integer N, minimize J2m=i ^m: where Si, S2,...,Sn 
range over all distinct A'"-tuples in A. 

Note that, if an element appears several times in A, it is allowed to appear the same 
number of times in the sum as well. 

Let us demonstrate the equivalence of the two problems. Given the sum X^^q^'"^* ' '^i; 
we may use the equality Wi • — Wi + Swi + 5wi -|- . . . -|- {2ni — l)wi to see that any feasible 
value for the objective function of the first problem is a feasible value for the objective 
function of the second problem as well. On the other hand, solving Problem 2 is trivial. 
Namely, one minimizes the sum there simply by taking the N least elements of the set A. 
In particular, for each i, the multiples of Wi present in the optimal solution will be all odd 
multiples Wi, Swi, 5wi, ... up to some {2ni — l)wi. Thus, the optimal solution of Problem 2 
yields the optimal solution of Problem 1 also. We note in passing that this discussion shows 
also that the minimum (for each of the problems) is obtained at a unique point unless A 
contains multiple elements. (However, we shall always refer to the optimal solution, even 
when there may be several.) 

A simple consequence of the above is 

Proposition 1. Let n = (rii)^^ be the optimal solution of Problem 1 for some value of N. 
Then the optimal solution of Problem 1, with N + 1 instead of N, is n' — (nQ^g^j where 
n'j = Uj + 1 for some < j < K — 1 and = rii for i ^ j. 

Remark. It is convenient to comment here on the effect a certain change in the original 
problem would make. One may consider the energies Ei to be WiUiini + 1) instead of Wiuf. 
This would change A to be the set of all even multiples of the Wj's. Obviously, this would 
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leave intact the equivalence of Problems 1 and 2. One can check that this would have also 
no effect on Theorems 1 and 2 infra. 

To formulate our main result we need a few definitions and notations. Real numbers 
91,02, ■■■ ,0r are independent over Q if, considered as vectors in the vector space R over 
the field Q, they are linearly independent. Equivalently, this is the case if the equality 
rriiOi + 171202 + . . . + rrirOr — for integer mi, m2, . . . , rrir implies mi = m2 = . . . = m^ = 0. 
Considering the actual physical system (a collection of metallic grains), it is reasonable 
to assume that the capacitors Cj are random, so that generically they are independent 
over Q. Without loss of generality we may rearrange the K capacitors such that Co = 
maxo<i<x-i Ci- It is also useful to divide all the capacitors by the largest one, so that the 
scaled capacitors Cj = Cj/Co with 1 = cq > Ci, C2 . . . , ck-i are dimensionless. Finally, set 

S = Co + Ci + . . . + Ck-1- 

Now we formulate our main results. 

Theorem 1. Suppose Co,Ci, . . . ,Ck-i are independent over Q. Then the sequence 
{Xn)n=i asymptotically F -distributed, where the distribution F is given by either of 
the following two representations: 

X <0, 

F{x) Q<x< 2wo, (4) 




2^0 < X, 



0, X < 0, 

E (1-51 + 1) n 11(1 0<x<2^o, (5) 
* s<z{i,...,K-i} ies i^s ^ "^"^0^ 

, 1, 2wo < X. 

It is not immediately obvious from the formulas, but F has one discontinuity, namely at 
the point 2wo- The reason is that, as the elements of A are all odd multiples of the Wj's, 
and as wq is the smallest of the WiS, it happens occasionally that there is no odd multiple 
of Wi, . . . ,wx-i between two consecutive multiples of Wq. The size of the atom at 2wo is 
1 ■ 11^1:1^(1 ~ Ci). This is easily explained intuitively. In fact, the "density" of odd multiples 
of Wi is Cj times the same density for multiples of Wq. Hence the "probability" that an 
interval of the form [(2n — 1)wq, (2n+ l)u'o) does not contain an odd multiple of Wj is 1 — Cj. 
Assuming that the "events" of containing different WiS are independent, we conclude that 
the proportion of multiples of wq in A whose successors are also such is n^i^(l ^ q)- Since 
the proportion of multiples of wq in A is 1 , we arrive at the required expression for the size 
of the atom. 
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Now we would like to study the asymptotic of the distances between consecutive elements 
of A as the number of capacitors grows. Obviously, as this happens, the distances become 
smaller. More precisely, on the average we have ^ odd multiples of each Wj in each unit 
interval, and hence we have there Y,f=o^ 2«r ~ 2w^ elements of A altogether. Hence the 
average distance between consecutive elements is To understand the asymptotics of the 
gaps, it makes sense therefore to normalize them so as to have mean 1. Thus, we multiply 
the distances by 2^;^, and ask about the asymptotic behavior 

Theorem 2. Suppose the capacitances Cq,Ci, . . . are chosen uniformly and independently 
in [0, 1]. For each K, let Fx denote the distribution corresponding to the normalized gaps 
when taking into account the hrst K capacitors only. Then, with probability 1, the distri- 
butions Fk converge to an Exp(l) distribution function. 

Remark. As will be seen in the proof, we actually use much less to prove Theorem 2 than 
is required by the conditions of the theorem. Namely, we need the capacitances Cj to be 
linearly independent over Q, and that they do not form a fast diminishing sequence. 

It is worthwhile mentioning that this type of "Poissonian" asymptotic behavior of con- 
secutive gaps is typical. For example, this is the case for uniformly selected numbers in 
[0, 1], and is conjectured to be the case in other interesting cases as well. (See, for example, 
pn| and |12[ and the references there.) 



In the course of the proof, we shall make use of the notion of uniform distribution 
modulo 1 and a few basic results relating to it. (The reader is referred to Kuipers and 



Niederreiter |]IU[ for more information.) A sequence {xnj'^^i of real numbers is uniformly 
distributed modulo 1 if 

|{l<n<M:a<K}<6}| 

— — > b — a, < a < b < 1 , 

where {t} is the fractional part of a real number t. In terms of Definition 1, (x„) is uniformly 
distributed modulo 1 if and only if the sequence ({a;„}) of fractional parts is F-distributed, 
where F is the distribution function of the uniform distribution on [0, 1]: 

{0, X < 0, 
X, < X < 1, 
1, X>1. 

The generalization of the notion of an asymptotically F-distributed sequence to that of 
an asymptotically well F-distributed sequence clearly carries over to our case. Instead of 
requiring only that the dispersion of large initial pieces of the sequence becomes more and 
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more even, we require this to happen at arbitrary locations. This version is termed well- 
distribution. Thus, {xn)'^=i is well-distributed modulo 1 if 

\{L<n<M:a<{xr,}<b}\ 

— > b — a, < a < b < 1 . 

M — L M-L^oo 

Both notions have multi-dimensional analogue. A sequence (x„)5^^ in R** is uniformly 
distributed modulo 1 in H.^ if 

\{l<n< N ■.a< {xn} <b}\ -A-/, x ^ . . 

~ ~ .r~ -^]lib^-a^), 0<a<b<l, 

1=1 

where inequalities between vectors in R'^ are to be understood component-wise, = 
(0, 0, . . . , 0) G R'^, a = (oi, 02, . . . , Os), and so forth. 

Perhaps the most basic example of a sequence which is uniformly distributed modulo 1 
is {na)^^i, where a is an arbitrary irrational. In the multi-dimensional case, the sequence 
{nai,na2, ■ ■ ■ ,nas) is uniformly distributed modulo 1 in if and only if the numbers 
1, ai, 0^2, . . . , are linearly independent over Q. Moreover, in this case uniform distribution 
implies well-distribution (cf. [rU, Example 1.6.1, Exercise 1.6.14]). 



Given a partition N = Uj=i and positive integers rj, j = 1, ...,/, we define the (rj)^.^j^- 
inflation of the given partition as the partition of N obtained by inflating each element of 
each of the sets Aj into elements. More precisely, we construct sets Bj, j = 1, ...,/, as 
follows. For a positive integer i, let f{i) = j if z G Aj. Given any positive integer n, let m 
be defined by EIILl^ /(^) < rn < f{i) . Let n E Bj if m E Aj. The following lemma is 
routine. 

Lemma 2. In this setup: 



1. If D(Aj) = dj, l<j<l, then D{Bj) = p'^' . 

2. IfBD{Aj) = d^, l<j<l, then BD{Bj) = [''^' . 



Proof of Theorem 1. Between any two consecutive odd multiples of Wq, there is at most 
one odd multiple of each wj, 1 < j < K — 1. In fact, one easily verifies that, given a positive 
integer m, there is an odd multiple of Wj between (2m — l)wo and (2m + l)wo, namely there 
exists an integer n with 

(2m - l)«;o < (2n - l)wj < (2m + l)wo, (6) 

if and only if 

l-Cj 1+Cj- 



mcj G 



(modi). (7) 
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Moreover, the relative position of {2n — l)wj within the interval [(2m — l)wo, (2m + l)wo) is 
the same, but in the opposite direction, as that of mcj (mod 1) within the interval (^-^, ^-^], 
that is 

(2n - l)wj = a • (2m - l)wo + (!-«)■ (2m + l)wo, (0 < a < 1) , (8) 
if and only if 

1 — c- 1 + c- 

mcj = (1 — a) • — ^-^ + a ■ — (mod 1) . (9) 

Next we define a partition of N as follows. Write the elements of A in ascending order: 
A = {6i < 62 < 63 < . . .}. Given n e N, let S" C {1, 2, . . . , - 1} denote the set of all 
those j's such that the unique interval of the form [(2m — l)wo, (2m + 1)wq) containing 6n 
contains an odd multiple of Wj. The set of all integers n giving rise in this way to any set 
S is denoted by Bs- Consider the partition N = [jsc{i,2,...,K-i} Bs- To prove the theorem 
using Lemma 1, we have to find the Banach densities of the sets Bs and the asymptotic 
distribution of the corresponding subsequences {Xn)n&Bs of Xn- 

The partition of N into sets of the form Bs is obtained as an infiation of a somewhat 
more straightforward partition. In fact, let S be any subset of {1,2, . . . , K — 1}. Denote 
by As the set of those positive integers n for which the interval [(2n — l)wo, {2n + l)wo) 
contains odd multiples of Wj for j G S* and does not contain such multiples of the other Wj^s. 
Then N = [jsc{i,...,K-i} is a partition, and its (15*1 + l)5c{i,2,...,_ft'-i}-iiiflation yields the 
partition N = U5c{i,2,...,i^-i} Bs- 

In view of the equivalence of (^) and (^), As is the set of those n's for which ncj G 
{-Y^, — ^ for j G S* and nCj ^ {~y^^ for 3 ^ S. By the conditions of the theorem, 
the numbers l,Ci, . . . ,ck-i are linearly independent over Q, and hence the sequence c = 
(nci, nc2, . . . , ncK-i)'^=i is well- distributed modulo 1 in R^~^. This means that 

D(A5) = BD(A5) = n 11(1 - c^) • (10) 
Denote the right hand side of (|IOD by ps- In view of the above and Lemma 2, this implies 

The denominator on the right hand side can be given a simpler form. In fact, let Xi, i = 
1, 2, . . . , ii" — 1, be independent random variables with ~ B(l, q), and X = Z^ill^ X^. 
Then: 

J2 {\T\ + 1)PT = E{X + l) = l + c, + ... + ck-1 = s . (12) 

TC{1,2,...,K-1} 
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Hence: 

BDiBs) = M±I)^ . (13) 

Let S be an arbitrary fixed subset of {1,2, . . . , K — 1}, say 5" = {1,2,...,/}, wliere 
< I < K — 1. IfnG As, tlien tliere exist odd integers ain, a2n, • • • , a^n sucli tliat ajnWj G 
[(2n - l)wo, {2n + l)wo). Put: 

Vn = {ainWi, a2nW2, • • • , ttinWi) - {2n - 1)1^0 ■ (1, 1, . . . , 1) G [0, 2woy, ne As . 

By tlie equivalence of (|) and (^, the sequence (fn)neAs is well-distributed modulo 2wo in 
Now each f„ gives rise to / + 1 terms of (Xn)neBg5 as follows. Let v^^^ < v^^ < • • • < v^p be 
all coordinates of f„ in ascending order. Set: 



The sequence (Xn)n€Bg consists of all coordinates of all vectors Un- Now we use the 
fact that if Xi,X2,...,Xr are independent random variables, distributed U(0, /i), and 
are the corresponding order statistics, then each of the random vari- 
ables - . . . , X('-) - X^^-^\h - has the distribution function defined by 
G{x) = 1 — {x/hy ioT < X < h (which follows as a special case from p|, p. 42, ex.23]). 
Consequently, for each 1 < j < / + 1, the sequence given by the jth coordinate of all 
vectors u„, n G As, is asymptotically well Gi-distributed, where Gi{x) = 1 — {x/2wq)^ for 
< X < 2wq. Hence the sequence (Xn)neBs is asymptotically well Gi-distributed. Combined 
with ([T3|) , it proves (H). 

We shall indicate only briefly the proof of (^, which is quite simpler. This time, we split 
(Xn) into a union of subsequences (Xn^); 0<i<-ft' — 1, by putting Xn in the sequence 
Xn'' if is a multiple of Wi. Clearly, the proportion of terms of (xn) belonging to (Xn^) is 
Ci/s. Next, consider the minimal odd multiples of all Wj's which are larger than The 
minimum of these K numbers is For each j ^ i, the distance from 5„ to the minimal 

odd multiple of Wj following 5„ is "distributed" U(0, 2wj). (For z = it is also possible that 
the next term will be again a multiple of wq.) The linear independence of the Cj's over Q 
implies that these K — 1 distances are (statistically) independent, so that their minimum is 
distributed according to the function G2{x) = 1 — 11^^^ (l ~ ^) on the interval [0,2wo). 
These considerations can be formalized to prove (^. This completes the proof. 

Remark. It is possible to shorten the proof by proving directly the equality of the right 
hand sides of @) and (|^). In fact, it is easy to integrate both forms with respect to x; 
the equality of the resulting expressions follows easily from the binomial theorem. We have 
chosen the long way, as it is more instructive. 
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Proof of Theorem 2. The distribution Fk is obtained from that in Theorem 1 by stretch- 
ing by the constant factor Hence: 



Fk{x) 



0, 



K-l K-1 



i=0 3=0 



X < 0, 



< X < s, 



s < X 



(14) 



Note that some of the values appearing on the right hand side depend on K imphcitly. 
Namely, since Wq is assumed in Theorem 1 to be the least Wj, each time a Ci is selected 
which is larger than all the heretofore selected Cj's, we have to rearrange the Cj's, thus 
changing and the Cj's. We have to show that 



FAx] 



K^oo 



l-e 



X > . 



Indeed, fix x > 0. Since 



s = Co + Ci + . . . + Ck-I 



Co + Ci + . . . + C, 



Co 



> Co + Ci + . . . + Ck-1 



and the Cj's are independent and uniformly distributed in [0, 1], we have 



a.s. 

S > OO . 

K^oo 



Hence, with probability 1, for sufficiently large K we have 

1 K-l K-l / 



FM = 1 



1=0 0=0 



X 

2^. 



Thus, to prove ([15|) we need to show that 



1 



K-l K-l 



1=0 j=o 



a.s, ^ ^ 



S J -ftT— >oo 

Now, on the one hand, using the inequality 



X > 



(15) 



(16) 



(17) 



(18) 



(19) 



l-t<e 



-t 



t G R , 



we have 



K-l 



J=0 



i = o,i,...,ir- 1 



and therefore 



K-l K-l 



i=0 3=0 



K-l 



i=0 



^^-x+x/s _ ^-x+x/s a.s.^ ^- 



(20) 
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On the other hand, as t — we have 



2 ' ' 2 

so that for all t in some sufficiently small neighborhood of 

e'(*+*') < 1 - t . 

Consequently: 



s 



2 



I 1 — ) > e J>^' > g ^ _ ^21) 



j=0 



Obviously, with probability 1, s grows linearly with namely for all sufficiently large K 
we have s > aK for a suitably chosen a > 0. (In fact, any a < \ will do.) By (pT]): 



j=0 



Jl ( 1 _ l^j > e— ^-^/^^ ^e-^ . (22) 



From dl) and (H) it follows that 



which completes the proof. 

A plot of P{x) is given in figure |l| and shows indeed that it has a maximum at x = 
and a delta function component at the inverse of the largest capacitance. In fact, when 
K oo, the weight of the delta function shrinks to zero and P{x) approaches the Poisson 
distribution e~^. 



II. FLUCTUATION OF COULOMB BLOCKADE PEAK SPACINGS IN LARGE 

SEMICONDUCTOR QUANTUM DOTS 

A. Motivation 

Recently, it became apparent that the physics exposed in the addition spectra of quan- 
tum dots is rather rich, and hence its investigation is at the focus of both experimental and 
theoretical studies. The present section concentrates on the distribution of spacings between 
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Coulomb blockade peaks in large semiconductor quantum dots. Coulomb blockade is evi- 
dently one of the hallmarks of mesoscopic physics. The experimental achievement of tracing 
an addition of a single electron to a quantum dot and the appearance of isolated conductance 
peaks led to the concept of single electron transistors. After the origin of Coulomb blockade 
peaks has been elucidated, investigation is directed toward more subtle questions like their 
heights, widths and spacings. The underlying physics is related to the ground state energy, 
chemical potential and inverse compressibility of a few electron island coupled capacitively 
to its environment, as well as fluctuations of these quantities with the number N of electrons 
on the dot. 

As far as the distribution of spacings between adjacent Coulomb blockade peaks is concerned, 
the question can be stated as follows: According to the simplest picture in which the quantum 
dot is regarded as a single electron island whose coupling with the leads is through its 
capacitance C, the total potential energy of a quantum dot is Q'^/2C — VgQ where Vg is the 
corresponding gate voltage. The conductance peaks occur at those values of Vg for which 
CVg = e{N+l/2) where e is the electron charge (henceforth e = —1). For this value of the 
addition of an electron to the dot (which contains electrons) does not cost any charging 
energy. The position of the A^*'^ Coulomb blockade peak is then a linear function of A^ and 
therefore, the spacing should be a constant 1/C, independent of A^. Recent experiments 
PJT3[] indicate however that spacing between Coulomb blockade peaks in small quantum 
dots is in general not constant but, rather, a fluctuating quantity close to Gaussian. The 
average of its distribution approximately coincides with the constant value mentioned above, 
but the elucidation of its standard deviation is still under investigation [0,|T^]. 
The situation is even less clear if the quantum dot is very large. As indicated in a series of 
recent experiments, the spacing occasionally vanishes, namely, two peaks (and sometimes 
even three peaks) coincide. Moreover, the evolution of peak positions and spacings with an 



applied magnetic field indicates the existence of strong correlations between them |[16|-[18 



These observation motivated numerous theoretical models based on the concept of pair tun- 
neling [|l^] or that of two-electron bound-states in depleted electron islands pO| . 
In the present section we examine the scenario according to which a large quantum dot like 



the one used in the last experiment ||T^, is in fact, composed of several electron islands 
which are coupled capacitively among themselves as well as to the leads. (This is a natu- 
ral extension of the scheme discussed in the first section in which no coupling is present). 
Electrons are added in such a way that the total potential energy of the dot is minimum. 
This simple generalization of the single island picture leads to a remarkable change in the 
spacing distribution from a Gaussian P,p^ centered around a finite average to a one which 
is large at small spacings. When the coupling between islands is weak, the distribution has 
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indeed a maximum at zero spacing. This result is short of explaining the perfect overlap 
of peaks, since it requires a delta function component at zero spacing. Yet, it leads to the 
occurrence of couples and sometimes triples of closely spaced peaks, similar to the experi- 
mental observation. Moreover, the evolution of the peak positions with the magnetic field 
is qualitatively similar to the experimental one. On the other hand, the present model does 
not predict a definite periodicity in the bunching of Coulomb blockade peaks with electron 
number N. In the next subsection the model is explained and the results of calculations are 
presented in the third subsection. 



III. FORMALISM 

Consider a large isolated two dimensional quantum dot in a perpendicular magnetic field 
B subject to a gate voltage Vg. Unlike the traditional Coulomb blockade picture it might 
contain several electron islands which can be regarded as metallic objects with inductive 
couplings among themselves. These are determined by a positive definite symmetric matrix 
C whose diagonal elements Cu = Cj > are the corresponding capacities whereas the non- 
diagonal elements Cy = Cji < , (z 7^ j) are the corresponding coefficients of induction. 
The electrostatic energy of such a system can be written as, 

Ec = lf: P^,N,N, - VgN, (24) 

where Ni is the number of electrons on island i (the number of islands K > 1 might be around 
10), = is the total numbers of electrons, and the (symmetric positive-definite) 

matrix p = C^^. 

Beside the electrostatic energy it is assumed that electrons in each island occupy single 
particle quantum states (orbitals) whose energies {i = l,2,..K;a = 1,2,..) depend on 
the confining potential as well as on the magnetic field. The latter is manifested through 
its orbital effects as well as due to Zeeman splitting (in which case the quantum number a 
contains also a spin label). The corresponding occupation numbers be either or 1. 

The system described above might then be represented by a classical Hamiltonian 

H = H, + H,p, (25) 



where the charging Hamiltonian He is just the electrostatic energy ^ written in terms of 
the orbital occupation numbers, 

K K 

Hc=^Y. Pij E ^ic] E -VgY.Y.^i», (26) 

^ i,j=l a a' i=l " 
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and the single particle part of the Hamiltonian, Hgp is, 

K 

The precise form of the matrix elements pij = [C^^]ij as well as the single particle energies 
eia are specified in the next section when we present our results. Despite the fact that the 
Hamiltonian H = Hc + Hgp is classical (and relatively simple) the elucidation of its spectrum 
for large and K is virtually hopeless. In order to compute the ground state energy E{N) 
one has to find the minimum of H on all the possible sets (A^i, N2...Nk) with the constraint 
Yl,f=i Ni = N. Note that the so called "Coulomb Glass" model obtains as a special case when 
i refers to a lattice site with random energy = ^i, and a single orbital Ni = rti = 0, 1. 
The interaction matrix is then given by pij = l/vij for i ^ j and pa = where rjj is the 
distance between sites i and j. 

The position of the A^*'' conductance peak is given by the first difference of the ground state 
energy, namely, the chemical potential of the (isolated) dot, /i = E{N + 1) — E{N). The 
spacing between peaks is determined by the second difference (the inverse compressibility), 
defined already in equation |I| 

The occurrence of close peaks for certain values of electron number correspond to 
small values of xn (recall that for a single island quantum dot in which the single particle 
energies are neglected, the inverse compressibility is a constant l/Cn). For the more general 
model described above the spacing distribution will of course fluctuate. In general, some 
constants appearing in the Hamiltonian H = + H^p are random, (e.g. the elements of 
the matrix C and the single particle energies eia), but most experiments are performed on 
a single quantum dot, so that fluctuations are meant with respect to the electron number 
A^. The numbers x{^) might then be considered as values assumed by a random variable 
X which has a certain distribution function P{x)- 



IV. RESULTS 

We now return to the full Hamiltonian of Eq. ^ It contains the coefficients of capaci- 
tance and induction matrix Cij and the single particle energies eia as input. In choosing the 
actual numerical values we use a few guidelines, one of them is to avoid too many indepen- 
dent input data. As will be clear below, these are not fitting parameters but rather, a set 
of constants which are chosen once for all on general physical grounds. First, for the elec- 
trostatic part, recall that the matrix C should be a symmetric positive definite matrix with 
Cii > and Cij < for i 7^ j. Besides, we expect it to be random. We then 
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are random numbers uniformly distributed between and W whereas the non-diagonal ele- 
ments are uniformly distributed between —w and 0. Evidently, in choosing these constants, 
one has to keep w « W in order to maintain positive definiteness. Actually, it is only the 
ratio w/W which matters, so one can assume W — 1 and use e^/W = 1 as an energy unit, 
leaving w as a constant reflecting the strength of coupling between islands. Second, for the 
single particle energies, we consider each electron island i as a two-dimensional potential well 
Vi(r) — \Mujfr'^^ where M is an effective mass. One may then regard OJ^^^"^ as a measure 
of the radius of the corresponding electron island. Since the capacitance Ca is also propor- 
tional to this radius we assume uji — ^C^"^, where 7 is a constant reflecting the relation 
between charging energies and single particle energies. The single particle energies in each 
electron island (subject to a perpendicular magnetic fleld B) are then known analytically. 
To be more speciflc recall that in two dimensions there are two quantum numbers for the 
orbital motion (denoted hereafter as n,m) to which we add a spin index a — ±1. Then, 
with a — {n, m, a), we have 



where Uc = eB/Mc is the cyclotron frequency, g is the g-factor and /is is the Bohr magneton. 
Since g contains the effective mass it is not known accurately. Its value is constrained 
on physical grounds (sec below). Note that in this scheme, the spacings between single 
particle energies arc deterministic, and do not follow the Wigncr surmise. The main cause 
of ffuctuation is then due to the combination of non-random single particle energies and 
the occurrence of numerous charging energies. Finally, the number of electron islands K 
is determined by the size of the quantum dot. The four input data of the model arc then 
K, w, 7 and g. Note that the gate voltage Vg does not have an important role here. Indeed, in 
actual experiments the variation of gate voltage serves to adjust the energies E{N) with the 
chemical potential of the leads but here the ground state energies are calculated directly. 
In order to avoid redundancy we are content with having a single set of these constants 
which is physically reasonable. In particular, it assures that the charging energy is much 
larger than single particle level spacings and that the Zeeman splitting is small at moderate 
magnetic fields. Specifically, we take w = 0.03, 7 = 0.1 and g is chosen such that for 
moderate fields the Zeeman splitting is of the order of the mean level spacing. 
With these prescriptions, we can find the ground state energy as a function of the magnetic 
field B for each electron number N. We use just a brute force trial algorithm, and hence 
cannot treat systems with large number of electron islands. An intriguing question is whether 
the addition of an electron leads to a redistribution of all electrons among the islands or else, 
the newly added electron will chose a place such that the energy cost is minimal while all 




(28) 
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others remain intact. Our results indicate that in most cases there is no redistribution, but 
in some cases redistribution does occur although it does not involve an overall re-shuffling. 
Consider as a representative example a quantum dot with K = b electron islands. The 
ground state energies E{N) and the inverse compressibilities x{^) (i^i the absence of a 
magnetic field) are calculated for electron number up to 200. The first question we 
addressed is how the electrons are added among the islands. Figure ^ shows electron 
numbers iVj in each one of the five islands, v.s the total electron number N, in the range 
60 < < 80. Evidently, the order of curves is according to the value of the capacitance 
Cii. On a larger scale, the numbers Ni grow linearly with N as it should be. Let us then 
consider the question of redistribution. From figure ^ we see that redistribution occurs 
only once, as grows between 67 and 68 (see the vertical line). The addition of electron to 
the second (or third) island involves also a transfer of an electron from the fourth island to 
the second (or the third) one. This scenario occurs also in other ranges of with the same 
proportion (namely about four percents). Thus, within the present model, redistribution is 
present although it is rare and minimal. 

The next question is related to the fiuctuation of the inverse compressibility x{^)- Figure 
^ displays x{^) v.s N for 40 < A^ < 80. From this figure one notices that the inverse 
compressibility fiuctuates quite strongly. In particular, it becomes small, (compared with 
its average) although it does not vanish. Note also that there is, at first glance, no traces 
of any periodic structure. This is also verified on applying its Fourier transform. The 
distribution P{x) is drawn in figure ^. It is indeed very different from the constant 
P(x) = S{x ~ l/C*) appropriate for a single island quantum dot with capacitive coupling C. 
A trace of this latter characteristic is the peak near x = 1. The eventual decrease of the 
distribution near x = seems to be related to the fact that, unlike the case of independent 
systems, the capacitance matrix is non-diagonal. Switching on the coupling here then has 
an effect similar to a weak "spacing repulsion" , similar to the familiar effect of perturbation 
in two level systems. Recall, however, that beside the capacitance induction matrix C the 
total energy is determined also by the single particle energies in each island. The remarkable 
point is that the distribution is not Gaussian (which is the hallmark of spacing distributions 
in small quantum dots). 

Finally, we check the behavior of the first difference E{N + 1) — E{N) as function of the 
magnetic field. Recall that this quantity is proportional to the position of the A^*'* Coulomb 
blockade peak. As a measure of the strength of the magnetic field we use the parameter 
uJc/ujq where Uq is the harmonic oscillator frequency of the largest island. The positions of 
the peaks for 39 < A^ < 48 are displayed in figure |^. Comparison with the results displayed 
in figure 2 of Ref. reveals a large degree of qualitative agreement. In particular, it 
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shows that groups of two (and sometimes even three) electrons can tunnel through the 
quantum dot at almost the same gate voltage. The oscillations at small magnetic field just 
mark transitions to lower Landau levels as the magnetic field increases. The phenomena of 
alternate bunching {N, N + 1) ^ {N — 1, N) is also reproduced in the present picture. 
In conclusion, we suggest a classical model in which a large semiconductor quantum dot 
is viewed as a collection of metallic electron islands with capacitive and inductive coupling 
among them. The effect of magnetic field is manifested through its orbital as well as its spin 
effects. The model can explain the occasional occurrence of couples or even triples of closely 
spaced Coulomb blockade peaks, as well as the qualitative behavior of peak positions with 
the magnetic field. 
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FIG. 1. Distribution of inverese compressibility P{x) for a system containing independent sub- 
systems each of which has a ground state energy proportional to the square of the number of 
particles it contains. The graph corresponds to a special case of equation ?? with K = 5. 

FIG. 2. Distribution of electrons among the five islands as function of the total electron number 
A'". The constants K, w, 7 and g are as explained in the text. The number Ni of electrons on 
island i is commensurate with its capacitance Ca. In most cases the addition of an electron does 
not perturb the occupation of other islands. An example of redistribution is marked with a vertical 
line. An addition of an electron causes a minimal redistribution. 

FIG. 3. Inverse compressibility xn as function of electron number in the dot. The units on 
the ordinate are the charging energies of the latgest island, namely, e^/W = 1 (see text). 

FIG. 4. Distribution P{x) of level spacings in the dot at zero magnetic field (not normalized). 
Note the peak at small values of %. 

FIG. 5. Coulomb blockade peak position for electron number N between 39 and 48 as function 
of the magnetic field. Here uJq is the oscillator frequency of the largest island and lOc is the 
cyclotron frequency. The units on the ordinate are energy units as explained previously. They are 
proportional to the gate voltage appropriate for the corresponding peak. 
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